Unified quaternionic description of charge and spin transport and intrinsic 

non-linearity of spin currents 
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We present a unified theory of charge and spin transport using quaternionic formalism. It is 
shown that both charge current and spin currents can be combined together to form a quaternionic 
current. The scalar and vector part of quaternionic currents correspond to charge and spin current 
respectively. We formulate a unitarity condition on the scattering matrix for quaternionic current 
conservation. It is shown that in presence of spin flip interactions a weaker quaternionic unitarity 
condition implying charge flux conservation but spin flux non conservation is valid. Using this unified 
theory we find that spin currents are intrinsically non linear. Its implication for recent experimental 
observation of spin generation far away from the boundaries are discussed. 
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During the past decade the study of Spin transport 
has emerged as an important area of research in con- 
densed matter. One of the important questions in spin- 
tronics is generation and manipulation of spin currents in 
semiconductors by means of electric field and spin-orbit 
coupling. Experimental observation of spin-Hall effect, 
both optically [H 0] and electronically 3 1 , acted 
alyst for recent theoretical works 0-|6|. Theoretically, 
the spin transport has been addressed using a variety 
of methods, e.g., drift-diffusion 0, QjKubo formalism (8|, 
and the scattering theory of transport (9rfl^|. The scat- 
tering theoretical formulation is particular suited for the 
mesoscopic system since it takes into account both quan- 
tum coherence effects and effects of measurement (i.e. 
of leads) into account. A detailed account of scattering 
theory for mesoscopic charge transport has been avail- 
able in Ref.[l3j]. The essential ingredients of scattering 
theory are formulated in terms of scattering matrix [l3|. 
For charge transport the charge conservation (global U(l) 
gauge symmetry) constraints demands that scattering 
matrix be unitary, i.e., SS^ = S^S = I. 

The situation for spin transport is rather subtle (see 
Ref . [lij ) . In presence of spin flip interactions (Spin-Orbit 
interaction, magnetic impurities etc.) though it is gener- 
ally stated that spin currents are not conserved however 
a mathematical equation describing spin flux non con- 
servation in terms of scattering matrix is still lacking. 
In the absence of such a guiding equation it is not clear 
which specific combination of the elements of scattering 
matrix should occur in the calculation of spin currents as 
opposed to charge current where unitarity of scattering 
matrix comes to rescue. This has lead to multitude of 
equation for spin currents in the literature (see Ref.0- 

On the other hand it is well know that scattering ma- 
trix of a SO coupled system is quaternionic^^, [la | re- 
flecting the break down of spin rotational symmetry. A 
quaternion q is defined as q = q° <E> 1% + iq ■ <x with 
(q°iQ) € C, where q is a vector and <t is a vector of Pauli 
matrices. The essential features of quaternion is that it 



combines the scalar and vector quantities into a single 
entity. The quaternion q can be viewed as having scalar 
and vector part(q = [Scalar, Vector]). Intuitively it is 
also clear that the charge currents are scalar in nature 
while spin currents are of vectorial nature. Therefore, 
it would be natural to look for a unified description of 
charge and spin currents in terms of quaternionic cur- 
rents. In the present article we show that both charge 
and spin currents can be combined together in terms of a 
quaternionic current, X — [I q , I s ], where scalar part gives 
usual particle(charge) current while the vector part cor- 
responds to the spin current. We further show that this 
naturally introduces a quaternionic metric on the Hilbert 
space. We formulate a Strong Unitarity condition for the 
quaternionic scattering matrix implying conservation of 
quaternionic metric and show that it is equivalent to con- 
servation of both charge and spin currents. Further it is 
shown that in presence of a spin flip interaction a weaker 
form of quaternionic unitarity condition is valid-implying 
particle flux conservation but spin flux non conservation. 
Using this weaker unitarity condition we derive equation 
for spin currents in two terminal system and study its 
implications. 

Before we proceed further it will be helpful to look at 
charge and spin currents separately. To elucidate the es- 
sential physics we confine our discussion to the case of two 
terminals systems. The Spin orbit coupled sample is of 
finite extent and is connected to reservoirs via the ideal 
leads. We assume that the leads have negligible spin- 
orbit coupling and spin relaxation so that spin currents 
are well defined in the leads far away from the sample- 
lead boundary. To describe the spin resolves channels in 
the lead we can choose a global spin quantization axis in 
all the leads along z direction ( spin resolved wavefunc- 
tions are eigen function of a%). The incoming and out- 
going spin resolved current amplitudes in leads is given 
by column vector A and B where A=[Al, Ar\ t (here T 
means transpose) and similar defining equation holds for 
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annihilation operator in channel m a) is a column vector 
with N^ a components corresponding to number of spin 
up and down channels. The B and A are related to each 
other via K x K scattering matrix (K = 2Nl + 2Nr) 

B = SA (1) 

Charge current or Particle Flux: The incoming and 
outgoing charge currents are given by the the bilinear 
combinations A^A and B^B. Mathematically these bi- 
linear combinations are in-fact inner product with respect 
to Identity metric, A given by, 

A= [{I 2 ®I L ,0},{0,I 2 ®/ K }] ee [{A £ ,0},{0,A fl }] (2) 

where J 2 being two by two identity matrix in spin space 
and II, Ir are Nz X Nl and Nr x Nr identity matrix 
for the left and right lead respectively(corresponding to 
spatial channel index). In terms of A, total incoming and 
outgoing particle (charge) fluxes can be written as, 

R 

II = Il,L+Il,R = A-^A=Y,AiA a A a 

a=L 

R.N a 

= ^ t ( a an a an Jra an a an) (^) 
a=L n=l 

R 

Ilut = IUl + IUr = B^AB = £ Bt A a B a 

R,N a 

= V (b u b a +tf- a b- a ) (4) 

a=L n=l 

Using Eq.(p} into Eq.(j4|) and comparing it with Eq.(j3]), 
we immediately see that charge flux conservation implies 
conservation of current metric A, i.e., 

S'AS = A. (5) 

In other words charge flux conservation conserves the 
identity matrix. In mathematical language it implies that 
on the Hilbert space the identity metric is preserved. 

Spin Currents or Spin Flux: In view of the above dis- 
cussion a natural question would be what is the corre- 
sponding metric for the spin currents? In other words 
how to combine A' ', A and B* , B to obtain incoming 
and outgoing spin fluxes respectively. Toward this end 
we define, the global spin current metric $7 whose matrix 
form is given by, 

n = [{a®I L ,0},{0,cr®I R }] ee [{n L ,0},{0,n L }} (6) 

where <x = {a x x + o~ y y + o~ z z) is vector of Pauli matri- 
ces, and I a is (N a x N a , N a being number of spin de- 
generate channels ) identity matrix in lead a. In terms 
of the spin current metric, incoming and outgoing spin 
fluxes are defined as, 



R 

i? n = A^nA = Atn a A a (7) 

a=L 

ts i is Ijs,x js,y xs,z-i /n\ 

± in,L ' ± in,R L J m s in ' in J V / 

R,N a 

= y ] [( fl an a a n + a a n a a n)i *( a <L n a an ~ a a n a un)i 
a=L n— 1 

( a an a an~ a an a an)] (9) 
R 

r out = B^fiB = Bin a B a (io) 

a=L 

— ± out,L i ± out,R — Y 1 outi 1 out^ 1 out\ I \ lL ) 

Where the explicit expression for outgoing spin flux can 
be obtained from Eq.([5]) by replacing incoming oper- 
ator a by corresponding outgoing operators b. From 
Eq.([5]) we notice that longitudinal spin current(spin cur- 
rent along the spin quantization axis) I?£ corresponds 
to the net spin polarization flowing in this is the stan- 
dard definition of spin current used in the recent works 0- 
H^ |. The spin currents 1^ and 1?^ are transverse spin 
currents (transverse to spin quantization axis) and corre- 
sponds to spin flip processes and deep inside the lead (far 
away from the sample-lead boundary) where SO coupling 
is absent such terms would vanish. However near the 
sample lead boundaries these will be non zero and reflects 
the quantum interference phenomena. These spin flip 
terms are similar to the complex mixing conductance of 
circuit theory developed by Brataas, Nazarov and Bauer 
[T3|. Using Eq.([T]) into Eq. ([T(If and comparing it with 
Eq. (JTJ), we at once see that spin flux conservation implies 
that the spin current Metric $7 is preserved, i.e., 

n = s^ns. (12) 

The above spin flux conservation equation actually a 
generalized vector unitarity condition. This is physi- 
cally plausible also since conservation of a vector quantity 
would imply both a conservation of its direction as well 
its magnitude. 

Quaternionic Currents: Now we show that charge cur- 
rent and spin currents can be expressed by a single equa- 
tion using quternionic formalism. Towards this end we 
define quternionic current metric in lead a is defined as, 
Q a = [A a ,n a ] = [J 2 <g> I a ,cr 55 I a ], and correspond- 
ing quaternionic current as X a — [2^ ; X s a ] . The quater- 
nionic current X a and quternionic current metric Q a 
can be viewed as having a scalar part and vector part 
corresponding to particle (charge) flux and spin flux re- 
spectively. The full quaternionic current metric for two 
leads system is expressed in matrix form as, 

Q = diag[Q L , ~Ql, Qr, -Qr) (13) 

By defining a combined vector creation operators as, 
P\ — {A^, B L , A r , B r ), we can write the net 
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quternionic current X = [X q , X s ], flowing through the sys- 
tem as, 

X = P^QP = J2[AlQ a A a -BiQ a B a ] (14) 

a 

= Y, A Uh ® la, n a ]A a - B\ t [I 2 ® I a , Sl a ]B&5) 

a 

It is easy to see that Eq. (fT5)l is nothing but a combined 
form of Eq.([TT]) and Eq.©. The Eq. (|T^j) is similar to 
the Eq.(l.ll) derived by Biittiker in Rcf. 13] for particle 
fluxes. We have show that it is valid for quaternionic 
currents as well. Now we discusses flux conservation 
for quaternionic current X. Using the defining relation 
Eq. ([1] in Eq. tjl 5|) we arrive at strong form of quaternionic 
unitarity condition, 

Q = S t QS=[A,n] (16) 

The above equation is nothing but a unified form of 
Eq.© and Eq. (fT2"| corresponding to charge and spin 
current conservation. It is straight forward to see that 
in presence of spin flip interaction the Eq.(|16p can not 
be satisfied for the vector part which would imply spin 
flux conservation as given by Eq. (fl~2|) . Towards this end 
we note that an element of scattering matrix can be ex- 
pressed as a quaternionic 

SafS = S% ® 12 + S*0 <7 X ~\~ B*p ® CT y + ® C7 Z (17) 

Since the spin flux conservation equationf (|12p implies 
that Scattering matrix commutes with the spin current 
matrix ft. This implies that each element of scattering 
matrix as given by Eq. (1171) commutes with Pauli matrices 
which obviously is not true. Therefore a weaker unitarity 
condition implying charge flux conservation but spin flux 
non conservation will be valid, i.e., 

S^QS=[A,X] = [A,S^nS]. (18) 

The Eq[l4]together with Eq[l8]provides a unified descrip- 
tion of charge and spin transport in presence of spin flip 
interactions. Using these two equation we find following 
explicit form of charge and spin current operators, 

£ = J, I dEY J Al{E)M M {a,E)A 1 {E) 

A-aSa/S^a^f Sl^p(E)A a Sa^ (E). (19) 



C = ^z\ dEYAl(E)r^(a,E)A,(E) (20) 

187 

T/3 7 = ttaifiapfiay ~ S a p(E)il a iS aj (E). (21) 

Where Q a % = o~i ® I a is the spin current metric along i 
spatial direction, where i = x,yorz. Note that the charge 



current operator (Eq. (fl"9)) ) is same as derived by Biittiker 
inRef.dl. 

Average Spin Currents: To find the average spin cur- 
rent we note that thermal averages of creation cb£ m and 
annihilation operator a G am for reservoir of unpolarized 
electrons is, 



4°m( E Kn( E ') = 6 (E - E')8 Pl 8 mn 5 aa ,fp (22) 

where fp = fp(E,fip) is the Fermi distribution function 
with the chemical potential \ip of reservoir /?[13|. Per- 
forming the average we find that the average longitudinal 
spin current I sz (along the spin quantization axis) for two 
terminal systems is, 

tsz " f J77i \ I jjcra i po" — a — T} — & — c\f 

a ~ An \\ H aa + H aa ~ K aa ~ H aa )Ja 

+ (Tap +T^° -T-°-°).M (23) 

with a = L, R and /3 ^ a. In writing down the Eq. (|23[) 
we have not assumed time reversal symmetry. In case 
time reversal is a symmetry, R a a a a = Raa~ a an d = 

T Pa~ a ^ T ap" = T Pa" stiU the E q-<E3|) predicts a non 
zero spin current. The important point to note is that 
even if one insists on using unitarity condition of particle 
flux sector (conservation of A in Eq. (fl8]l ) still the spin 
currents can not be written down as a difference of Fermi 
distribution functions. To see this explicitly we rewrite 
the eq. ([23|) using conservation of A as, 

la = J dE ^(Ka" - C)/» + TaflUp ~ fa)] , 

(24) 

where T a p = (T£| + T°~° - T~™ - T^~") and f t = 
fi(E,Hi) for i — a or j3. From eq.fJM]) we see that 
net longitudinal spin currents has equilibrium as well 
nonequilibrium contribution and these can not be sep- 
arated. Because for the equilibrium part the spin re- 
solved reflection probabilities contribute while for the 
non-equilibrium part corresponding transmission proba- 
bilities contribute. Experimentally one measures the net 
currents (net spin polarization) and not the transmission 
or reflection probabilities separately. Therefore spin cur- 
rents are intrinsically non-linear in an electrical circuit 
where the driving force is electrical voltage. This is not 
surprising because the gauge conjugate to spin currents 
is non-Abelian SU(2) gauge which consists of terms de- 
scribing spin-orbit interaction see Ref.dlGJ]. While the 
Eq. (j23|) is written in terms of U(l) gauge (electrical volt- 
age) conjugate to charge currents. We stress that the 
non-linearity of spin currents does not imply breakdown 
of linear response theory, rather it is non-applicability of 
linear response with respect to electrical voltage for spin 
currents. 

Transverse Spin Current: The expression for average 
transverse spin currents can not be written down in terms 
of spin resolved probabilities. It turns out to be much 
more complicated, however, to bring out essential physics 
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we express these currents in terms of spin resolved block 
scattering matrix as, 

(25) 

In the above equation a and ft can take values either L 
or R corresponding to left and right lead. Similarly the 
summation over a' can take two values a or — a. It is 
obvious from the Eq. ([25|) that it can not be expressed 
as a difference of Fermi distribution function signifying 
the non-linearity as discussed above. In Eq.([23|) the spin 
conserved and spin flip scattering matrix elements (reflec- 
tion and transmission amplitudes) occur together which 
is a quantum interference effect. On thermal averaging 
such an effect would vanish deep inside the lead far away 
from the sample lead boundary. However near the bound- 
ary this will give rise to oscillatory spin currents where 
evanescent state plays a role [ill EH- As remarked ear- 
lier these currents are similar to the mixing conductance 
which appear in the circuit theory of Brataas, Nazarov 
and Bauer (l7j . 

As we have seen above that longitudinal spin currents 
survives deep inside the lead and may give rise to spin 
currents far away from the boundary. In a recent ex- 
periments Sih. et al.Q observed spin generation away 
from the edges of a GaAs sample subjected to an elec- 
trical field. A theoretical explanation of such an effect 



has been worked out in Ref.0-@- In Ref.Q it was ar- 
gued that spin generation far away from the boundary 
is not possible in the linear regime. Though we can not 
make a direct comparison with the experimental results 
of [2|, however, the Eq.(l23| for longitudinal spin current 
above generically show that these are intrinsically non 
linear with respect to electric field hence can exist far 
away from the boundary. 

Finally let us discusses the possibility of observing 
these experimentally. Most promising technique looks 
to be an optical measurement, especially using Magneto- 
optical Kerr effect(MOKE). In the Kerr effect the ro- 
tation of the polarization axis of reflected beam is 
proportional to the spin polarization along the beam 
direction [2(j. Therefore to study all the component of 
spin current it would be essential to measure Kerr ro- 
tation with different angle of incidence (beam direction) . 
Considering that such experiments are routinely done for 
spin studies it should not be very difficult to experimen- 
tally test the effects discussed in this paper. 

To Conclude we have presented a unified formulation 
of spin and charge transport using quaternionic currents. 
Unitarity condition of quaternionic currents is formulated 
and spin currents equation are derived. It would be inter- 
esting to extend this formulation to study spin shot noise 
which might have some nontrivial bearings on semicon- 
ductor spintronics. 
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